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Abstract 

By superposition of regular gauge instantons or merons, ensembles of gauge fields are constructed 
which describe the confining phase of SU(2) Yang- Mills theory. Various properties of the Wilson 
^ loops, the gluon condensate and the topological susceptibility are found to be in qualitative agree- 

^ ment with phenomenology or results of lattice calculations. Limitations in the application to the 

xy^ glueball spectrum and small size Wilson loops are discussed. 
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Introduction 

Despite convincing evidence for confinement in Yang-Mills theories from lattice gauge 
calculations and despite a host of analytical studies and conjectures, the intricacies of non- 
abelian gauge fields have prevented an unambiguous identification of the mechanism of 
confinement. In particular, the gauge dependence of many of the mechanisms proposed 
constitutes a severe obstacle in reaching this goal. I will choose a diflFerent path towards 
understanding confinement. I will describe construction and analysis of ensembles of gauge 
fields which exhibit confinement and account for other facets of the dynamics of SU(2) Yang 
Mills theory. The studies to be presented are based on the ideas that instantons [Ij or 
merons [2j, solutions of the classical Yang-Mills equations, play a prominent role also in the 
quantum theory i4j . They are related to investigations which led to the instanton gas and 
liquid models O E] which successfully describe a variety of phenomena of strong interaction 
physics [71 E]. These models however miss confinement, the trademark of the dynamics 
of Yang-Mills fields. I will discuss an alternative construction of ensembles of gauge-fields 
obtained by superposition of "regular gauge" instantons or of merons [9j . Although irrelevant 
for single instanton properties, the choice of either the regular or the singular gauge in the 
superposition of instantons leads to two different phases, the gas or liquid phase of singular 
gauge instantons with a well defined low density limit and the strongly correlated confining 
phase of regular gauge instantons (or merons). 



Gauge fields of instantons and merons 

After appropriate choice of the coordinate system in color space and after regularization 
of the singularity, the SU(2) gauge field for a meron or an (regular gauge) instanton in 
Lorenz gauge with the center at the origin, is given by 

with the Pauli matrices a^, the 't Hooft tensor rja/^u [7J, and with <^ = 1 denoting merons and 
C = 2 instantons. The size parameter p controls the short distance behavior of instanton 
and meron fields. For instantons, and for merons with vanishing p, a^(x) is a solution of the 
Euclidean classical field equations [K . 

Important for the following are the infrared properties of the constituents. The gauge 
fields of a single meron and of a single regular gauge instanton behave asymptotically as 1 /x 
while the asymptotics of the field strengths, 

F;;AA]=d,A:-d^A; + e'^''^AlAt, 

is different for merons and instantons 

X ^ oc , F^^i ^ , F^=2 , (2) 

due to the cancellation of the abelian and non-abelian contributions to leading order for 
instantons. Therefore the action density 

s{x) = UvF,,F^^ (3) 
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behaves asymptotically as x ^ and x ^ for merons and instantons respectively and gives rise 
to an infinite action for merons. The topological charge density is defined by 

~s{x) = \ivF,,F^\ (4) 

and yields a finite value for the topological charge 

p = i/2. (5) 

The fast asymptotic decrease of the field strength of an instanton is the basis for the alternate 
representation of instantons in "singular gauge" 

sing Vafiiy^iy /^x 

where the gauge field decays asymptotically as and f] denotes the anti-self-dual 't Hooft 
tensor. 



Gauge fields from superposition of merons or instantons 

The SU(2) gauge fields to be considered are superpositions of instanton or meron fields 

Np 

A^{x) = - z{i))h-\i). (7) 

1=1 

The dynamical variables of these fields are the positions of the centers z{i) of the instantons 
or merons and their color orientations 

h{i) = h^{i) + ih(z) . cr , hl{i) + = 1. (8) 

In general, the action of gauge fields generated by superposition of either merons or in- 
stantons is logarithmically divergent in the infrared. By a judicious choice of the color 
orientation, finite values of the total action can be obtained. This is achieved, e.g., by the 
following superposition of 4 instanton or merons 

3 

A{x) = ^ Rla{x - z{i)) + a{x - z{A)) , 

with the color rotations RJ by tt around the three color directions i. With this choice 

^ ^ oc : ^ ^ , F[a] - ^ , s{x) - ^ , (9) 

the system is "neutral" and the action is finite. The constituents, instantons or merons, are 
"confined". Removal of a single constituent gives rise to an infinite action. Thus ensem- 
bles of instantons and merons exhibit similar infrared properties and in both cases strong 
correlations between the constituents are required to guarantee a finite action. 
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The ensembles of gauge fields (cf. Eq.([7|)) to be constructed in the following are defined 
by the partition function 

Z = / dz,d/i,e"^^'^^'^'^^^' , (10) 



with the effective action S identified with the standard (Euclidean) action (cf. Eq. ([3])) 

S = [ d^xs{x). (11) 
Jv 

This definition of the partition function guarantees that for finite a non-zero weight is 
assigned to neutral configurations only, i.e. to fields with the asymptotic behavior Q. After 
construction of ensembles of field configurations using the Metropolis algorithm, vacuum 
expectation values of observables O can be computed 

(O) = I y dzidhie-^''^^^'-'''^^0[A{zi, hi)] . (12) 

In the numerical determination of these ensembles, the location z{i) is restricted to a hy- 
percube 

-l<z,ii)<l, (13) 

an equal number of instanton and anti-instantons (merons and antimerons) is used, and in 
most of the applications the values of size and coupling constant are chosen as 

p = 0.16, ^ = 32. (14) 

Characteristics of instanton and meron ensembles 
Action densities 

In this section I will characterize qualitative properties of gauge fields in the ensembles 
defined by (10). In Fig. [l] is shown a typical landscape of the action density displaying peaks 



of single merons on top of a background. Background and single meron contributions to 
the total action are of the same order of magnitude. Unlike in the superposition of singular 
instantons, the appearance of single meron or instanton structures is not automatic, it rather 



refiects the dynamics defined by the partition function (10). This becomes evident in the 
comparison (Fig.jl]) with the action density of a meron field of the corresponding stochastic 
ensemble {g'^ = oc). The values of the action density of the two ensembles 

g^ = 32, {s) = 1500, 200 < s{x) < 11000, 

/ = (X) , (s) = 110000 , 15000 < s{x) < 330000 , (15) 

differ by almost two orders of magnitude. As a result of of the huge background generated 
by the stochastic superposition of the constituents no single meron contribution (s^ax = 
18000) can be identified. The maximum close to the center of Figjl] is a result of the 
stochastic superposition as can be verified in an analytical calculation. The reduction of the 
action from the strong coupling value is a first quantitative signature of the importance of 
correlations in the {g'^ oo) ensembles. Another measure of the strength of the correlations 
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FIG. 2: Response of the action of 50 instantons to changes of the color orientation along three 
independent geodesies on the manifold of the color rotations ([s]) (left) of the instanton closest 
to the center and to changes of its position (right). 



between constituents is the response of the action to changes of the color orientation or 
of the position of a single meron or instanton \L0\. As shown in Fig.[2| under variations 
of the color orientation of a single constituent the total action is increased by up to 60% 
and by up to 20% under variations of the position. In a weakly interacting system of 50 
instantons (instanton gas) one would expect a 2% decrease of the action when removing a 
single instanton and even smaller changes under variations of the color orientation. 



5 



Confinement 



In this section I will present evidence that the phase of highly correlated constituents is 
confining. In Yang Mills theories, the area law of Wilson loops is the signature of confine- 
ment. The Wilson loop is defined by 

W = ^trjPexpi^ ^ A^{x)dx^^ , (16) 

with the integral ordered along the closed path C. The following results refer to rectangular 
paths with aspect ratio 2:1. For sufficiently large loops the results for the logarithm of the 
expectation value of the Wilson loop can be parametrized as a sum of constant, perimeter 
(V) and area (A) terms, 

ln{W) =uj + tV-(jA, (17) 

with a denoting the string tension. The parameters of the fit for ensembles of 500 merons 
(p = 0.16) are 

UJ = -2.1, T = 2.3, a = 22.9 (18) 
for the stochastic ensemble {g'^ = oc). For ensembles with g'^ = 32 the following values 

cj = -0.7, T = 1.1, a = 11.8 (19) 

are obtained if the color orientations are dynamical variables while the positions are ran- 
domly chosen but fixed, and 

cj = -0.7, r = 1.1, a = 11.5 (20) 

for the ensemble with dynamical color orientations and dynamical positions. These three 
ensembles give rise to an area law with values of the string tension differing by a factor 
of 2 between stochastic and dynamically correlated ensembles. The large difference in the 
dimensionless ratios formed by action density and string tension, 

^^ 500,(^2 = 00), ^^^n{g' = 32), (21) 

reffects the different dynamics of the ensembles. As in lattice gauge theories, confinement in 
the strong coupling limit is a result of disorder of the system generated by the unconstrained 
fiuct nations of the gauge fields. Suppression of the fiuctuations if g'^ ^ oo affects more 
strongly the local action density than the non local- Wilson loop and in turn leads to the 
strong decrease in the ratio g'^{s)/a. Only minor changes of Wilson loops are obtained 
if, in addition to the color orientations, also the meron positions are treated as dynamical 
variables. This result is in accordance with the above findings concerning the response to 
changes in the color orientation and position (Fig. [2]) of single constituents. It suggests that 
the confining phase is a close relative of the nematic phase of liquid crystals with strong 
correlations in internal and comparatively weak correlations in position space pT] . 



Quantitative results 
Scaling and confinement 

Quantitative predictions for observables evaluated in the ensembles defined above require 
definition of a scale. As in lattice gauge calculations the phenomenological value, a = 
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4.4 fm ^, of the string tension is chosen to set the scale. Rescahng the different contributions 



to the Wilson loop (17), 



ln{W) =uj + rVxV -aXA, (22) 

and fitting the scaling parameter A separately for each ensemble gives rise to a universal 
curve as shown in Fig.|4} The values of the universal parameters and of the unit of length 
are for the instanton ensembles 

Lu = -0.52, T = 1.18, a = 19.0, 1 u.l. = 2.08\/A fm . (23) 

Similar results have been obtained for the meron ensembles. 

As the Wilson loops, also other observables exhibit approximate independence of the number 
of degrees of freedom. Tablejl] summarizes for ensembles with Nj instantons the results for 
the vacuum expectation value of the action density s and the topological susceptibility (also 
in physical units), 

(si)'/ ^'^(^»^~(0)) ' (24) 

which is a measure of the strength of the fiuctuations of the topological charge In 
comparison to the values obtained from sum rules ^3 and from lattice gauge calculations 

M 

85 fm"^ < {s) < 260 fm"^ , (25) 

the expectation value of the action density i.e. the gluon condensate is of the correct order 
of magnitude. It receives about equal contributions from the background field generated by 
the superposition of the constituents and from the peaks of the fields of single constituents 
as suggested by Fig.[TJ Also the value of the topological susceptibility is in reasonable 
agreement with the lattice result x^^^ ^ 215 MeV [14j. The topological susceptibility is 
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FIG. 4: Logarithm of a Wilson loop as a function of its area for instanton ensembles Nj = 50 (star), 
100 (diamond), 200 (square), 500 (circle). The area has been rescaled with A given in Table [l| Also 
shown is the curve corresponding to the parametrization (22) with the values of the parameters 
given in (23). 
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TABLE L Properties of instanton ensembles: Vacuum expectation value of the action density s 
defined in ^ and the topological susceptibility x ( [24| ) (also in physical units), with the standard 
choice of the parameters ( 14) for ensembles of Ni instantons and with instanton density n/ = NijLr" 
and the values of the scaling parameter A (cf. ([22])) and the size parameter p in physical units. 



dominated by the topological charge density of single constituents. The single constituent 



contributions in the evaluation of the integral (24) yield the result, 



1/4 



0.83 ny^ 



1/4 



which agrees within 10% with the result of the numerical evaluation (TablejT]). Due to the 
dominance of single instanton properties similar values of this observable in the liquid phase 
[7J of singular gauge and in the confining phase of regular gauge instantons are obtained. 
Table[l] shows that ensembles with increasing numbers of constituents (Nj) and located in 
a hypercube of fixed volume (13) describe, after rescaling, systems of essentially the same 
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FIG. 5: Distribution of the values (cosi9) of Wilson loops for an ensemble of configurations con- 
taining 500 merons with sizes 0.48 x 0.24, 0.72 x 0.36, 1. x 0.5 in comparison with the distributions 



density rij of instantons or merons. The increase in number is accounted for by the increase 
in volume of the hypercube. The value of the string tension determines the density of 
constituents. Variations in the constituent size in the interval 0.06 < p < 0.25 have been 
considered and found to leave the scaling properties essentially unperturbed. Origin of the 
scaling properties is the scale independence of the asymptotics of regular gauge instantons 
or merons (Eq.[T]). In ensembles of singular gauge instantons on the other hand, p controls 
the strength of the asymptotics of these fields (Eq.[6]) and therefore the overlap and in turn 
the strength of the interaction of these constituents [16j . 

Wilson loops and glueball masses 

In addition to the area law, also more detailed properties of Wilson loops have been found 
to agree with those of lattice gauge calculations. I mention the Wilson loop distributions 
shown in Fig. [s] which, as similar lattice results [17j, can be interpreted as distributions of a 
diffusion process on S'^, the group manifold of SU(2), given by 



with the diffusion time t determined by the expectation value of the Wilson loop. Diffusion on 
entails "Casimir scaling" of the Wilson loop expectation values in higher representations 
which has been observed in lattice gauge calculations p^. Also in agreement with lattice 
gauge calculations is the breakdown of Casimir scaling observed for sufficiently large loops 
for ensembles of gauge fields similar to the ones discussed here . 

In a study of correlation functions of circular loops of equal radius r which are parallel to 
each other and orthogonal to the direction of separation (t) 



(26). 




(26) 



n=l 



Crit) = (H;(n,xo,t)H;(n,xo,0)) 



(27) 
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indications for formation of flux tubes in connection with the area law have been found. If 
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FIG. 6: Left: Logarithm of correlation functions of circular Wilson loops for an ensemble of 500 
instantons of size p = 0.08 as a function of separation in time t for 3 values of the Wilson loop 
radius r — 0.25 (black circles), r = 0.32 (red squares), r = 0.40 (blue diamonds) . Right: The 
same as a function of the variable rt. 



gauge strings are important degrees of freedom the application of the Wilson loop operator 
to the vacuum should generate such a gauge string with energy E ^ 27rcrr, and the relevant 
variable for describing the correlation function should be rt. Indeed the three correlation 
functions are described approximatively by a universal curve in terms of this variable as 
shown in Fig. [6} A rough estimate of the string tension based on this calculation yields 
a= (3.7±L3)fm"^ 

As the last topic I discuss the glueball spectrum calculated via (Euclidean) correlation 
functions 

{O{x)O{0)) - J2{n\0\n)e-^-''{n\0\n) . 

n 

of appropriately chosen operators carrying the quantum numbers of the glueballs. For 
sufliciently large separations, the correlation function is dominated by the state of lowest 
energy which can be excited by the operator O. For determination of the masses, local and 
non-local operators have been used and various ensembles differing in the size parameter and 
the number of constituents have been studied. The final result is shown in Fig.[7| Qualitative 
agreement with respect to the ordering of the level has been reached. The scale however 
is not compatible with the scale set by the string tension. The values of the masses have 
been multiplied with a factor of L8 (L5 for meron ensembles) to reproduce the 0+ glueball 
mass obtained in lattice gauge calculations. Determined by the requirement of the infrared 
finiteness of the action, the degrees of freedom do not adequately account for the dynamics 
on the scale of the size of the glueballs (0.2 -0.5 fm)[21j. Missing strength in the ultraviolet 
also shows up in the Wilson loops of small size as is shown in the right part of Fig.[7| 
While the stochastic ensemble with its unconstrained fiuctuations generates a Coulomb-like 
behavior at small distances, the ensembles of dynamical merons or instantons fail in the 
description of Wilson loops in the perturbative regime. A remedy for the shortcoming in the 
calculation of the glueball spectrum could be the inclusion of perturbative contributions in 
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FIG. 7: Left: Spectra of glueball states. Lattice gauge theory [20j (left part) Instanton ensembles 
after rescaling (right part). Right: Enlarged portion of Fig. [3] for small size Wilson loops. Upper 
curve obtained for the ensemble with dynamical meron positions and color orientations, lower curve 
for the stochastic ensemble {g^ — oc). 



addition to the non-perturbative correlation functions, as applied for ensembles of singular 
gauge instantons [22j. 

Conclusion 

The central achievement of these studies is the construction of ensembles of gauge fields 
which exhibit confinement. These fields are obtained by superposition of regular gauge 
instantons or merons. It is remarkable that superposition of the gauge equivalent regular 
and singular instantons gives rise to two distinctively diflFerent phases, the gas or liquid phase 
generated by singular and the confining phase by regular gauge instantons respectively. The 
confining phase owes its existence to the non-trivial requirement of finite action in the 
presence of the slow asymptotic 1/x decay of the regular gauge instanton or meron fields. 
Unlike the disordered strong coupling limit, the confining phase is an ordered phase similar 
to the nematic phase of liquid crystals. The infrared behavior of the constituents is not 
only the source of confinement it's scale independence also induces the scaling properties 
of observables under changes in the number of constituents. Besides various properties of 
Wilson loops, such as the area law, Casimir scaling or fiux tube formation, also the action 
density (gluon condensate) and the topological susceptibility are reasonably well reproduced. 
These latter quantities are determined by the relative strength of the gauge fields close to the 
center of the constituents and of the background gauge field generated by the superposition 
of the constituents. If observables such as the topological susceptibility are dominated by 
single instanton properties similar results in the liquid and in the confining phase can be 
expected. Limits in the applicability of meron and instanton ensembles are encountered if 
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short distance properties are important as is the case for the glueball spectrum and for small 
size Wilson loops. 
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